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We introduce a systematic low-energy approach to strongly correlated electron 
systems in infinite dimensions, and apply it to the problem of the correlation-induced 
metal-insulator transition in the half-filled Hubbard model. We determine the low- 
energy scaling functions of the metallic state, including the single-particle Green 
function and dynamical spin susceptibility, as well as thermodynamic properties and 
relate them to experimental data in transition metal oxides. 
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The correlation driven metal-insulator transition is a fundamental problem in condensed 
matter physics which has attracted renewed theoretical and experimental interest. From a 
theoretical point of view, the problem is fascinating because of its non-perturbative nature, 
the absence of a small parameter and the emergence of new low energy scales. Early ideas of 
Mott |TJ], Hubbard Brinkman and Rice [[5] have recently been put on a more quantitative 
footing by the development of a mean-field approach to the strong correlation problem which 
becomes exact in the limit of large lattice coordination The work of several groups |] [?J 
has established the existence of a Mott transition in the paramagnetic phase of the half-filled 
Hubbard model on frustrated lattices, which is driven by the collapse of an energy scale, 
the renormalized Fermi energy. The presence of this small energy scale made it difficult to 
analyze the critical behavior near the transition. 

In this letter, we introduce a novel projective self-consistent approach to solving corre- 
lated electron problems in large dimensions at low energies. We use the natural separation of 
energy scales to our advantage and are able to extract, for the first time, exact information 
about the critical behavior at the transition in infinite dimensions. We calculate the single 
particle spectral function near the Mott transition, which can in principle be measured in 
photoemission experiments. We establish that the coefficient of the u 2 term in the imagi- 
nary part of the self-energy diverges as the square of the coefficient of the linear term in the 
specific heat and, with few additional assumptions, relate it to the observed T 2 resistivity 
in the LaSrTiO^ system. We show that the local spin susceptibility diverges in the same 
way as the linear term in the specific heat and calculate their ratio, a generalized Wilson 
number. Finally, we show that the linear coefficient of the imaginary part of the dynamical 
spin susceptibility, relevant to both NMR and neutron scattering experiments, diverges as 
the square of the linear coefficient of the specific heat. 

Besides these results of theoretical and experimental interest, our method can be gen- 
eralized to deal with other problems, in which the separation of energy scales typical of 
strongly correlated electron systems can be exploited |]]. Two examples are the breakdown 
of Fermi liquid theory in strongly correlated electron systems |§ and superconductivity 
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in models with repulsive interactions [|10yiT|| , which have been recently analyzed using the 
large-d techniques. 

The Hamiltonian of the Hubbard model is given by 

h = - E (tn + tfdftfi* + £(/U - \)Ulh - ~) (i) 

<ij>o i Z Z 

where the chemical potential /x is taken to be zero, and the hopping matrix is scaled as 
tij — > t/y/d. In the limit of infinite dimensions, all local correlation functions of the lattice 
model can be calculated in terms of a local action, 

S imp [G ] = - f drdr'Y,fl{r)G-\r-r')Ur') 

+ [dTU{n n {r)- l -){n n {r)-\) (2) 



provided that G Q satisfies a self-consistency equation [O. In order to make contact with 



physical systems in finite dimensions, for which the bandwidth is finite, we consider the 



case of a semi-circular bare density of states, p {e) = (2/ttD)^1 — (e/D) 2 for which the 
self-consistency equation has the form 

G-\iuj n ) =iuj n + /j- {D/2) 2 G{iuj n ). (3) 



The local action (g) can be rewritten in terms of an Anderson impurity model p2 



Ham = Y, e k c l*Ck,* + Y, V k(flck* + h.c.) + U{n n - -)(n n - -). (4) 

ka kcr ~ ~ 

The self-consistency equation (|3|) is satisfied provided that the bath dispersion and the 
hybridization coupling Vk are determined through 

T,- ±L — = G^n). (5) 

where G{iuj n ) = — Jq dre %WnT < T T f a (r)fl(0) >h am ^ s the impurity Green function. 

For sufficiently strong interactions, a separation of energy scales occurs as demonstrated 
explicitly by previous numerical work on the Hubbard model ||5|-0. It has been shown that 
the one particle spectral function can be decomposed into a sum of a low and a high energy 



part, p(e) = p low (e) + p hl9h (e). p low (e) contains all states up to a cut-off that we take to be 
the Kondo temperature or renormalized Fermi energy of the half-filled Hubbard model and 
carries spectral weight A 

A = £V fe 2 /£ 2 (6) 

k 

where the primed summation runs over the low energy states only p h% 9 h [e) describes the 
upper and lower Hubbard bands, two atomic-like features at energy scales ±~, and carries 
spectral weight 1 — A. When A is small the presence of the low energy scale AD makes the 
system of equations (f|) and (||) numerically untractable. 

The main idea of the projective self-consistent approach is to eliminate the high energy 
states to obtain a low energy effective problem involving p low only and thus containing only 
one energy scale. To carry out this program we separate the impurity configurations of the 
impurity Hamiltonian (^) into a low energy sector f^\E >, with eigenenergy —U/4, and a 
high energy sector \E > and fjfj\E >, with eigenenergy U/4, where \E > denotes the empty 
configuration. Furthermore, the self-consistency equation allows us to divide the conduction 
electron bath into three bands as illustrated in Fig. ([!]): a metallic band centered around the 
Fermi energy, and semiconducing valence and conduction bands centered around energies 
±£7/2. For the analysis of low energy properties, the valence and conduction bands can be 
taken as dispersionless [13|, with the corresponding atomic states created by rf_ and 77+, 



respectively. 

The low energy effective Hamiltonian can be derived through a canonical transformation. 
In addition to the high energy impurity configurations, which can be eliminated by the 
standard Schrieffer- Wolff [ITJJ transformation, the self-consistency requires the elimination 



of the valence and conduction bands of the electron bath. For this purpose we first separate 
out the atomic Hamiltonian 
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Hatomic = U(n n - -){n n - -) 

+ ^(V+,aV+,a ~ V-,aV-,,) + ^ E + h.C.) (7) 



where V is determined from the fact that both the upper and lower Hubbard bands have a 
spectral weight (1 — A), which together with the self-consistency equation (El) implies that 



V = (D/2) yj (1 — A)/2. The lowest-energy eigenstates of this atomic Hamiltonian are a spin 
doublet \a > which correspond to the impurity spin doublet fl\E > dressed by the valence 
and conduction electrons. All other atomic configurations are located at energies higher by 



at least ~ ^. 



We have carried out the canonical transformation analytically to order (D/U) 2 ||T3| . |T5 
resulting in a low-energy effective Hamiltonian H e // = e s H.e~ s of the form 

He// = X! Jkk'S • Skk 1 + € k c ka C ka (8) 
kk' ka 

where J kk > = 8^^(1 + - |^A), S = \ YL aa > X^a^i and s kk > = \ J2aa' 4,<t^<«t ,c *V- 
Here, X a/3 = |a >< where \a > and \f3 > are the eigenstates of the atomic Hamiltonian, 
are projection operators. 

Heff corresponds to a Kondo problem coupling the atomic spin doublet to the low energy 
conduction electron bath. The low energy part of the Green function can now be calculated 
directly from H e ff, 

G lm »{iu n ) = - J dre^ < T T F(r)Ft(0) > He// (9) 
where F a = e s f a e~ s is the canonically transformed single particle operator and has the form 

F a = J2'a k [(Xaa ~ X^ a )Cka + 2X a ^C k -a] (10) 
k 

with ak = ^(1 + f fjr — ffjrA). The self-consistency equation for the low energy Green 
function then becomes 

k lU) n — €k 

The projective self-consistent method thus results in the closed set of equations (|B]-|TT|) which 
form the basis of our low energy analysis. The system contains only one energy scale, 
AD, which allows us to define rescaled variables, V k = Vfe/(\/A£)), e k = ek/(AD), iuo n = 
iu n /{AD), and H e// = E eff /(AD). 



In the following, we will focus on the critical properties of the Mott transition approached 
from the metallic side, for which it suffices to analyze Eqs. (|8|-|TT|) to leading order in A. 



Restricting to leading order in A implies that the requirement that the rescaled spectral 
weight is unity, / dCjp low {uj) = 1, can be satisfied together with the self-consistency equations 
(P|-|TT|) only for one particular value of U. This is precisely the critical interaction U c . 

The problem is now reduced to a self-consistent determination of the low-energy parame- 
ters ik and \4, which is carried out iteratively using a recently introduced zero-temperature 
algorithm |l6|Jl7[| . The functional equations (§) and flllD are truncated by representing 



the electron bath in terms of a finite set of free electron orbitals q, i = 1,...,N — 1, 
with the corresponding energy levels and hybridization matrix V*. The ground state 
of the impurity problem, now defined on an N site cluster, is obtained by the modified 
Lanczos technique f[8| |. The local Green function as well as correlation functions are cal- 



culated using a continued fraction expansion []T9| , |T6|] . The projected self-consistency con 



dition (prTJ) can be implemented in terms of a x fitting of the Matsubara Green function 
X 2 = Y2™=h min I E/c i^r k ~ < F] F >H cff (^n)l 2 where VL min and VL max are the low and high 
frequency cut-offs, respectively JT7|] . 

We have solved the equations for clusters of 4,6,8 and 10 sites. The critical interaction 
U c is found to be essentially unchanged as the number of sites is increased from 4 to 10, 
varying from 3.03 to 3.04 f20fl . In the inset of Fig. (|2|) we show the single particle spectral 
function for N = 10. The large number of poles clearly demonstrates the power of the 
projection onto the low energy sectors. The Matsubara Green functions as functions of the 
rescaled Matsubara frequency for N = 6, 8 and 10 are shown in Fig. (|2|). The low frequency 
part of the Green function improves systematically as N increases, and remains essentially 
unchanged for N = 8 and 10. Also shown in Fig. (Q) is the Green function for the non- 
interacting problem. The value of the scaled Green function at zero frequency is determined 
by the low energy self-consistency equation and is the same as that of the Green function of 
the non-interacting problem, consistent with the pinning of the density of states expected 
from Fermi liquid theory |2]J • The solution for the Green function also allows us to determine 
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the scaled self-energy T, low (ioj n ) by extracting from £(zu; n ) = iu n — (D/2) 2 G(iu n ) — G~ x (iu) v ) 
the terms with a singular dependence on A. Both G low (iuj n ) and 'E low (iu) n ) can be fitted by 
quadratic polynomials at low frequencies with the expressions 



G low (iCj n )= -i2.0 + 3.3(iw n ) + i3.9(iu n ) 2 } 

T} ow (iCj n )= -1.7iu n - il.l(iu) n ) 2 . (12) 

The term of the self-energy linear in iu n implies a quasi-particle residue z = (1 — 
c?£/ 'diuo n )~ l = A/ 1.7 which vanishes as the critical point is approached. The momentum- 
independence of the self-energy in turn leads to a quasi-particle mass m* /m — 1/z — 1.7/A, 
and a linear coefficient of the specific heat, 

7 ~ M (13) 

which diverge at the critical point. This divergence is consistent with the Brinkman-Rice 
scenario of the Mott transition 0], as well as the previous result in the Hubbard model 
at half-filling within second order perturbation theory . The term of the self-energy 

quadratic in (iu) n ) 2 gives rise to an imaginary part of the analytically-continued self-energy 

/m£(^ + *0 + ) = -H^ (14) 

which also diverges at the critical point. 

The local dynamical spin susceptibility Xsii^n) = (fi^s/^) 2 Jjf dTe luJnT < 
T T S z (r)S z (0) > Heff , where S z = |e 5 (/|/ T - /j/Je" 5 = (X TT - X u )/2, can also be cal- 
culated from a continued fraction expansion. The result in terms of rescaled frequency, 
Xs(i&n), is shown in Fig.(|3|). At low frequencies, Xs(i&n) can be fitted by 

X l r^n) = (W2) 2 (9 - 47p n |)/AD (15) 

The local static susceptibility is given by the constant term of Eq. (p~5|), 

Xs = (W2) 2 9/(AD) (16) 
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from which we derive a generalized Wilson ratio at the critical point, 

/ free 

R _ XslXs loc 2 (17) 

<-y /jfree y ' 

where x^ioc = 16/(37rD)( a ^ § -) 2 and 7/ re e = j^^b-i are the static local spin susceptibility and 
the linear coefficient of the specific heat for the free electron gas with a semicircular density 
of states, respectively. Compared to the universal value for the infinite bandwidth Anderson 
impurity model, Ram = 2, the critical value R is enhanced as a result of the finite bandwidth 
of the electron bath. In large dimensions the susceptibility at a generic q is not affected 
by the antiferromagnetic exchange interaction and, hence, has the typical value X\ oc - The 
enhancement of R is consistent with the ferromagnetic tendencies found in the Gutzwiller 



approximation |25j where the magnetic exchange is ignored as well. We note that the usual 
definition of the Wilson ratio in a lattice model is given in terms of the q = component of 
the static susceptibility, x s [f\ = 0). The antiferromagnetic exchange interaction does affect 
Xs(q = 0) (as well as Xs(q) for special q points in the Brillouin zone), leading to a vanishing 
ratio Xs(q = 0)/7 at the transition point. [^6| 

Finally, the term of the imaginary part of the dynamical spin susceptibility linear in 
\uj n \ implies that lim^ Xs (^-bo + ) = 47/(-DA) 2 (g/i#/2) 2 . If we use xioc for the static spin 
susceptibility this implies a finite generalized Korringa ratio at the critical point, which is 
again modified from the universal value for the infinite bandwidth Anderson model |27[ . 



Our results have direct implications for the Mott-Hubbard systems. The conductivity 
can be estimated by converting the uo 2 at zero temperature to (vrT) 2 at finite temperature 
and assuming that the self-energy we derived at the critical point also applies to the case 
of a hypercubic lattice in d dimensions [28[|. Using the Kubo formula, we find a resistivity 



p(T) = AT 2 where A = — ~575^^~pr giving rise to a finite ratio 

— = 8.2 x 10 3 a (rim) (18) 

where a is the hypercubic lattice spacing in units of meter. For Sri_ x La x TiOz with x = 0.9, 
which is close to the Mott transition, Eq. ( |18|) yields A/7 2 = 4.4 x 10 _6 f2m when we 
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use a = 8.3 x l(r 10 m This result is in good agreement with the measured value 

A/7 2 = 8 x 10~ 6 Qm p| . 



In summary, we have introduced a projective self-consistent approach to strongly cor- 
related electron systems, which allows us to determine the critical properties at the Mott 
transition in the half-filled Hubbard model. We find low energy scaling functions that can 
be quantitatively related to experimental results in transition metal oxides. 

We thank V. Dobrosavljevic, A.E. Ruckenstein, and X.Y. Zhang for useful discussions. 
This work was supported by the NSF under grant ^DMR 9224000. 
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FIGURES 

FIG. 1. Schematic plot of the spectral functions of the conduction electrons and the impurity 
configurations, illustrating the low and high energy parts. 

FIG. 2. Imaginary part of the Matsubara Green function versus the rescaled frequency Co n for 
system size N = 6, 8, 10 and non-interacting (semi-circular) density of states. The inset shows the 
single particle spectral function as function of the rescaled frequency to with a broadening 5 = .01 
for N = 10. 

FIG. 3. Real part of the Matsubara spin susceptibility as a function of the rescaled frequency 
Cj n for N = 10. 
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